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Past  research  on  the  Bomberger  lotsize  scheduling  problem  has  relied  heav- 
ily on  the  following  assumptions: 


1.  Each  product  is  to  be  produced  on  a regular  invariant  cycle. 


2.  The  inventory  level  at  the  onset  of  each  production  run  for  a product 
being  produced  is  zero. 
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Abstract 


A Linear  Model  for  the  Single  Machine 
Lotsize  Scheduling  Problem 

by 

Thom  J.  Hodgson 

Past  research  on  the  Bomberger  lotsize  scheduling  problem  has  relied 
heavily  on  the  following  assumptions: 

1.  Each  product  is  to  be  produced  on  a regular  invariant  cycle. 

2.  The  inventory  level  at  the  onset  of  each  production  run  for  a 
product  being  produced  is  zero. 

In  this  paper  these  two  restrictions  have  been  dropped.  It  is  shown  that 
the  inventory  cost  function  can  still  be  expressed  as  a linear  function.  A 
linear  programming  model  is  formulated.  Sensitivity  analysis  is  used  to 
find  an  optimal  cycle  length.  Numerical  examples  are  used  to  show  how  the 
model  can  be  used  to  develop  lotsize  schedules  which  are  both  practical  and 
low  cost. 
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= Set  up  cost  of  product  i ($) 
c^  = Cost  per  piece  of  product  i ($/unit) 
d^  = Demand  rate  of  product  i (pieces/unit  time) 

F(T)  = The  inventory  carrying  cost  per  unit  time  for  a particular 

lotsize  scheduling  problem  as  a function  of  the  cycle  length  T. 

h = Inventory  carrying  charge  ($/$inventory/unit  time) 


I. 

J 


= The  inventory  level  of  product 
run  j . 


at  the  beginning  of  production 


k. 

3 


= The  index  of  the  product  produced  on  the  jth  production  run  of 
the  schedule 


N = Number  of  products 

n^  = The  number  of  times  product  i is  produced  in  the  schedule 
p.  = Production  rate  of  product  i (pieces/unit  time) 
s^  = Setup  time  of  product  i (time) 

T = Length  of  the  total  schedule  (time) 


W.  = The  idle  time  immediately  after  the  completion  of  the  jth 
J production  run  of  the  cycle 

X.  = The  start  time  of  the  jth  production  run  of  the  cycle 

z.  = The  index  of  the  next  production  run  of  product  k.  after  pro- 
'*  duct  ion  run  j. 


Introduction 


The  lotsize  scheduling  problem  in  its  various  forms  is  a problem  that  has 
plagued  both  practitioner  and  researcher.  As  its  name  implies,  it  has  elements 
of  both  inventory  and  scheduling.  The  objective  function,  the  sum  of  setup  and 
inventory  costs,  is  derived  from  inventory  theory;  and  the  feasibility  require- 
ment associated  with  machine  availability  is  primarily  a scheduling  problem. 

The  problem  is  basically  that  of  finding  a minimum  cost  production  schedule  on 
a single  machine,  for  N different  products,  each  with  known  demand  rate,  pro- 
duction rate,  and  setup  time.  The  problem  can  be  approached  in  two  different 
ways.  It  can  be  approached  as  an  infinite  horizon  problem  where  the  objective 
is  to  determine  a cyclic  repeatable  schedule  which  minimizes  setup  plus  inven- 
tory costs  per  unit  time.  It  can  also  be  approached  as  a finite  horizon  pro- 
blem where  the  objective  is  to  determine  a schedule  over  some  planning  period 
which  minimizes  setup  plus  inventory  costs  over  that  period.  This  paper  is 
concerned  with  the  former  approach. 

Rogers  [15]  provided  an  early  comprehensive  treatment  of  the  problem.  He 
structured  the  basic  problem  and  treated  the  problem  of  schedule  feasibility  in 
great  detail.  Bomberger  [3]  developed  what  might  be  considered  the  first  algor- 
ithmic approach  (Dynamic  Programming)  for  the  problem.  In  addition  (and 
possibly  more  important),  he  developed  a standard  data  set  of  10  products 
which  has  provided  the  yardstick  by  which  most  of  the  following  researchers 
have  evaluated  their  efforts.  The  literature  contains  a succession  of  papers 
[1,  2,  4 - 16],  each  providing  some  additional  insight  into  the  structure  of 
the  problem  as  defined  by  Bomberger. 

Tiie  rather  unique  and  insightful  approach  of  Elmaghraby  and  Mallik  [5] 
deserves  comment.  They  structured  the  problem  using  a network  model.  The 
power  of  the  approach  lies  in  its  conceptual  simplicity  and  its  adaptability 
to  both  finite  and  infinite  horizon  problems.  An  efficient  computational 


1 


procedure  for  the  approach  has  yet  to  be  developed. 


It  should  be  noted  that  the  Bomberger  definition  of  the  problem  does  not 
necessarily  match  what,  in  fact,  happens  on  the  production  floor.  A basic 
assumption  is  that  each  product  is  produced  at  a given  constant  interval,  and 
that  the  inventory  for  each  product  reaches  zero  just  as  its  production  run 
begins.  As  a consequence,  the  sum  of  setup  and  inventory  cost  per  unit  time  as  a 
function  of  the  cycle  length  T for  a product  is  just  [10] 

~ + he  Cl  - d/p)d  (1) 


where  A = setup  cost  ($), 

h = inventory  carrying  charge  ($/$  inventory/unit  time), 
c = piece  cost  ($/unit) , 
d = demand  rate  (pieces/unit.  time),  and 
p = production  rate  (pieces/unit  time). 

The  least  cost  cycle  length  under  these  assumptions  is  just  [10] 
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hc(l  - d/p)d 


(2) 


liquations  (1)  and  (2)  are  widely  known  and  used.  However,  it  has  been  this 
author's  experience  that  the  basic  assumptions  associated  with  the  cost  structure 
of  equation  (1)  arenot  always  valid  when  one  works  on  the  production  floor.  Machine 
availability  constraints  may  force  the  use  of  unequal  Jotsizes,  staggered  cycle 
patterns,  and  nonzero  inventories  at  the  beginning  of  production  runs.  This 
phenomenon  was  recognized  by  Maxwell  [13,  14]  in  his  early  work.  Baker  [1] 
also  recognized  the  practicality  of  unequal  lot  sizes.  Maxwell  formulated 
the  lotsize  scheduling  problem  as  a mathematical  program.  His  model  requires 
a repeatable  sequence  of  jobs  as  input  with  each  product  (i)  being  produced 
n.  times  during  the  sequence.  The  formulation  has  linear  constraints. 
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However,  the  objective  is  quadratic  and  non-convex.  In  general,  the  mathe- 
matical program  is  not  readily  solveable. 

In  succeeding  sections  of  this  paper,  a model  requiring  a repeatible  se- 
quence of  jobs  as  input  is  formulated  in  which  staggered  cycle  patterns  and 
non-zero  inventories  at  the  beginning  of  production  runs  are  allowed.  The 
model  has  linear  constraints  and  a linear  objective.  Consequently,  it  can 
be  optimized  using  linear  programming.  It  is  shown  that  the  optimal  length 
of  the  schedule  (cycle  time)  can  be  found  using  sensitivity  analysis.  An 
example  is  given  to  show  how  the  model  can  be  used  to  find  useful  production 


Model  Formulation 


Consider  the  following  situation.  Product  i (with  demand  rate  d^)  is 
produced  on  a single  machine  n^  times  every  T time  units.  Each  time  it  is 
produced,  d^T/n.  units  are  produced.  It  is  desired  that  product  i be  setup 
and  produced  exactly  every  T/n^  time  units.  The  desired  (least  cost)  inventory 
function  is  shown  in  Figure  la.  However,  product  i must  compete  for  the 
available  time  on  the  machine  with  the  other  N - 1 products  assigned  to  the 
machine.  In  order  to  achieve  a feasible  schedule,  it  is  deemed  necessary  to 
move  the  schedule  for  a production  run  of  product  i ahead  in  time.  The  inven- 
tory function  for  this  situation  is  shown  in  Figure  lb.  As  can  be  seen  from 
Figure  lb,  the  inventory  level  only  reaches  the  level  I at  the  stage  of  the 
rescheduled  production  run  and  is  inflated  over  part  of  the  next  inventory  cycle. 
This  results  in  an  increased  inventory  cost  for  the  schedule  which  is  proportional 
to  the  shaded  area  in  Figure  lb.  Since  the  lotsize  schedule  will  repeat  itself 
every  T time  units,  the  increase  in  total  inventory  cost  per  unit  time  can  be 
calculated  by  evaluating  the  area  of  the  shaded  area  (Figure  lb),  multiplying 
by  the  inventory  carrying  cost  (h)  and  the  piece  cost  (c^),  and  dividing  by 
the  cycle  length  (T): 

A cost  = h • c.  • (shaded  area)/T. 

The  shaded  area  is  easily  seen  to  be  equal  to  E • G , where 
C = (1  - di/pi)d  .T/n^  and 
E = I/(d.(l  - d./p.)). 

Therefore,  multiplying  and  simplifying, 

A cost  = h • c,  • I/n^  (3) 
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Note  that  the  increase  in  cost  per  unit  time  is  linearly  dependent  on  the  minimum 
inventory  level  I. 

With  this  information  in  hand,  let  us  proceed  to  a linear  formulation  of 
t lie  N-produet,  single  machine,  lotsize  scheduling  problem.  Suppose  that  a 
sequence  K (of  length  L)  of  production  runs  to  be  performed  during  a cycle 


(of  length  T time  units)  is  given,  is  the  index  of  the  product  to  be 


produced  on  the  j th  production  run  of  the  cycle.  Let  X^(0  < < T)  be  the 


start  time  of  the  jth  production  run  of  the  cycle,  > 0 be  the  amount 


of  idle  time  immediately  after  the  completion  of  the  jth  production  run,  and 
n.  be  the  number  of  times  the  index  i appears  in  the  sequence  K.  Note  that 


any  production  run  for  product  i lasts  d^T/p^n^  time  units.  The  following 


set  of  equations  relate  the  start  times  (X.)  and  idle  times  (W.). 

3 3 


VT 


X + 


+ W.  + S,  = X. 


J VV  j kj+i  “j+1 


3=1,  ....  L - 1 


(4) 
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Let  z.  be  the  index  of  the  next  production  run  of  product  after  production 


j.  Note  that  during  any  production  run  for  product  i.  d^T/n^  units  are 


produced.  The  following  set  of  equations  define  the  inventory  levels  (1^.) 


for  each  product  at  the  beginning  of  each  production  run. 


f + dk  T/nk  - d (X  - X ) = Iz  3 = 1,  2,  ... 
J 3 3 3 3 3 


(5) 


Note  that  equat.^n  (5)  is  not  necessary  when  production  run  j is  the  nk  th 

3 


traduction  run  for  product  k.  since  only  n - 1 equations  are  necessary  to 

3 Rj 


fully  define  the  relationships  between  the  nk  inventory  levels. 


The  objective  of  the  linear  program  is  to  minimize  the  total  inventory 
cost  (setup  cost  will  be  dealt  with  at  a later  point).  It  should  be 


«_ 'l 


* ■'  * j * 


6 


W ' 


- — ■; — v 


not'd  that  if  all  of  the  inventory  levels  (1^)  are  zero,  the  only  inventory 


isLs  are  just  the  balanced  cycle  inventory  costs. 


h } c.(l  - d./p.)d,T/2n. 

/.I  ill  i 

i=l 


Anv  additional  inventory  costs  due  to  modifying  the  balanced  cycle  schedule 
to  achieve  feasibility  are  accounted  for  by  equation  (3).  The  objective 
function,  then,  is  as  follows: 

N L 

min  fh  [ V c. (1  - d./p.)d.T/2n.  + 7 c,  I./n  ]}  (7) 

. , i xii  l .S  k.  l k. 

i=l  J=1  3 J 

Combining  the  objective  (7)  with  the  constraints  (4)  and  (5),  and  including 
Liu  usual  non-negativity  constraints  constitutes  a linear  program. 

Two  points  can  be  made  immediately  concerning  the  characteristics  of 
optimal  so1  "ions  to  the  linear  program.  First,  it  is  clear  at  least  one  of 
the  minimum  inventory  levels  (1^)  f°r  each  product  (i)  will  be  equal  to  zero.  This 
is  easily  proven  by  contradiction  and  is  left  to  the  reader. 

Second,  if  T is  included  as  a decision  variable  in  the  linear  program 
and  no  bounds  are  placed  on  it,  the  optimal  solution  will  occur  when  T is  at 
i!  minimum  feasible  value.  The  minimum  feasible  value  of  T can  be  calculated 

as  fci  I I ows  : ,, 

N 

} n . s . 
x i 

T . = -1--- :. . (8) 

min  N 

1 -/.VPi 
1=1 

This  second  characteristic  is  difficult  to  prove,  because  as  T decreases 
tlx  minimum  inventory  levels  (I.)  increase  at  an  unspecified  rate.  However,  it 


Determi ning  the  Optimal  Value  of  T 


1!  the  feasibility  constraints  are  ignored,  the  least  cost  value  of  T can 

lie  determined  from  the  following  formula: 

1/2 


T*  = 


r n 

2 }’  n.A. 

i=l  1 1 

N 

h } c . d . ( 1 - d./p.)/n. 

‘•-ii  1*1  1 


i = l 


l l 


(9) 


II  the  feasibility  constraints  are  considered,  it  is  easy  to  see  that  the  optimal 

* 

value  of  T will  be  no  less  than  T . Therefore,  as  a first  step  to  determining 
iIm-  optimal  value  of  T,  the  following  constraint  can  be  added  to  the  linear 


pn  ram: 

T max{T*,  Tminl  (10) 

In  order  to  determine  the  optimal  value  of  T,  it  is  necessary  to  determine 
the  optimal  value  of  the  objective  function  of  the  linear  program  (inventory 
cost)  for  all  values  of  T.  As  it  turns  out,  this  can  be  done  in  a straight 
forward  fashion  using  sensitivity  analysis.  Once  the  inventory  carrying  cost 
function  F(T)  has  been  determined,  the  value  of  T which  minimizes  the  sum  of 
setup  and  inventory  costs  can  be  determined. 

Assume  that,  for  a particular  set  of  products  and  a given  sequence  of 
production  runs  to  be  performed  during  the  cycle,  the  linear  program  has  been 
formulated  (Including  constraint  (10) ) and  solved.  Assume  that  in  order  to  main- 
tain feasibility,  it  is  necessary  to  schedule  jobs  such  that  several  of  the 
starting  inventories  are  positive. 

T may  he  increased  up  to  a certain  value  and  the  optimal  basis  will  still 
he  maintained.  The  amount  that  T can  be  increased  with  no  change  in  basis  can 
he  determined  by  the  entries  in  the  column  of  the  surplus  variable  for  constraint 
(10)  (column  v)  of  the  final  tableau  (See  f 1 7 ] , page  129-133,  for  a lucid 


eyp  I anat  i «>n  ) . 


Let  I . be  the  basic  variable  in  the  linear  program  solution  associated  with 

row  r of  the  final  tableau.  Then  the  tableau  entry  (a  ) in  the  rth  row  of  the 

surplus  variable  column  represents  the  amount  1^  will  change  if  T changes  by  one 

unit  (since  the  surplus  variable  is  represented  in  the  initial  tableau  as  a 

negative  unit  vector,  the  sign  of  the  entry  (ary)  is  reversed).  By  taking  the 

ratio  of  the  inventory  level  to  the  tableau  entry  (I^/a^)  for  each  row,  the 

range  of  T,  within  which  the  optimal  basis  is  unchanged,  can  be  determined.  In 

addition,  the  sensitivity  of  the  objective  function  to  T is  contained  in  the 

objective  row  of  the  surplus  column.  As  a consequence  of  the  above  discussion,  - 

it  is  seen  that  by  solving  the  linear  program,  a segment  of  the  inventory  carrying 

cost  function  can  be  determined  (at  this  point  it  should  be  apparent  that  the 

inventory  carrying  cost  function  is  piecewise  linear),  along  with  bounds  on  the 

range  of  the  sement.  If  T is  increased  beyond  the  range  of  the  segment,  the 

optimal  basis  changes,  and  a new  segment  of  the  inventory  carrying  cost  function 

can  be  established.  In  general,  when  the  optimal  basis  changes,  one  or  more 

of  the  basic  I.  will  leave  the  basis,  and  a like  number  of  non-basic  W.  will 

j J 

enter  the  basis.  Since  the  I^'s  decrease  as  T increases,  removing  I^'s  from  the 
basis  and  replacing  them  with  VL's  which  are  not  included  in  the  objective  will 
increase  the  sensitivity  of  the  objective  function  from  segment  to  segment.  In 
other  words,  the  inventory  carrying  cost  function  F(T)  is  convex. 

The  problem  of  finding  the  optimal  value  of  T,  then,  is  just  the  problem 
ol  minimizing  the  sum  of  the  piecewise  linear,  convex  inventory  carrying  cost 


function,  F(T),  and  the  hyperbolic  setup  cost  function 


Thu  graph  in  Figure  2 serves  to  illustrate  the  problem.  The  minimum  total 
aist  will  occur  either  at  one  of  the  breakpoints  or  in  between  breakpoints.  To 


determine  if  the  minimum  occurs  between  two  breakpoints,  the  following  procedure 
can  he  used.  The  expression  for  the  variable  elements  of  the  total  cost  function 
in  a given  segment  is  just 

N 

Total  Cost  = y n.A./T  + ST,  (12) 

i=l  1 1 

where  S is  the  sensitivity  to  T of  the  linear  program  objective  function  in  the 

segment.  The  minimum  of  (12)  occurs  at 
r N I 1/2 

| y n . A . 
ii 

r J1-^ d3) 

L S 

if  T'  falls  within  the  boundaries  of  the  segment,  then  the  minimum  occurs  within 
the  segment.  A simple  procedure  can  now  be  given  for  finding  the  optimal  value 

for  T: 

1.  Set  the  lower  bound  on  T according  to  equation  (10)  and  solve  the  linear 
program.  Determine  the  bounds  on  T for  the  linear  segment. 

Solve  for  T'  using  equation  (13).  If  T'  is  greater  than  the  upper  bound 
for  the  segment  C.0  T0  3.  If  T*  is  within  the  bounds,  ST0P:  T'  is  optimal, 

if  T’  is  less  than  the  lower  bound,  ST0P.  The  lower  bound  is  optimal. 

1.  Increase  the  Lower  bound  on  T to  the  upper  bound  of  the  segment  and  update 
the  linear  program  to  its  new  optimal  solution.  Determine  bounds  on  T 
for  the  new  segment.  G0  T0  2. 
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An  Example 


Ear  an  example  let  us  turn  to  Bomberger's  [3]  famous  ten  product 
data  set  which  Is  found  In  Table  1.  Doll  and  Why  bark  [4]  have  pre- 
sented the  lowest  cost  solution  to  date  for  this  data  set  ($32 .071/day) . 

With  a total  schedule  length  (T)  of  187.395  days,  products  1 through  10  are 
produced  1,  4,  4,  8,  4,  2,  1,  8,  4,  and  4 times  respectively  during  the  cycle. 

A sequence  of  40  jobs  in  Table  2,  which  satisfies  the  above  cycle  frequencies, 
was  used  lo  generate  the  linear  program  model.  As  expected,  the  procedure 
described  in  the  previous  section  achieved  the  same  cost  as  Doll  and  Whybark 

k 

(since  the  sequence  is  feasible  at  T , no  sensitivity  analysis  was  required). 

!U>v . ver,  the  solution  is  of  limited  practical  value.  It  is  difficult  to  imagine 
nr.  company  that  would  plan  to  produce  products  every  23.424,  46.848,  93.696, 
or  187.395  days.  It  boggles  one's  mind  to  think  of  the  reaction  a production 
foreman  might  have  to  such  a schedule. 

in  fact,  it  is  not  an  uncommon  practice  [7]  for  companies  to  limit  production 
cycles  to  specific  lengths  (i.e.,  one  week,  two  weeks,  one  month,  etc.)  which 
can  he  easily  implemented  by  production  personnel.  In  the  present  case, 
ii  1 1 1 g the  value  of  T at  1 year  (240  days),  and  fixing  the  individual  cycle  lengths  of 
the  10  product,  example  at  1,  3,  6,  or  12  months  (assuming  a 20  working  day 
month)  should  result  in  a more  useable  schedule.  However,  it  is  impossible 
to  develop  a feasible  schedule  using  conventional  means.  By  shifting  production 
run.,,  tiK  linear  program  is  able  to  find  a feasible  solution  for  a cost  of 
‘•12. 365/d  iv  (see  Table  3 for  the  L.  P.  schedule).  For  less  than  a 1%  increase 
in  setup  ind  inventory  costs,  the  model  is  able  to  provide  a workable  schedule. 

• lie  point  is  rli.it  a production  planner  can  use  the  linear  programming  model  as  a 
tool  for  t to  development  of  good  practical  schedules.  The  model  does  optimally 
what  production  planners  and  schedulers  can  only  now  do  in  a heuristic  fashion. 


Table  1 

Boinberger's  Data* 

.Setup 

Setup 

Piece 

Production 

Demand 

Time 

Cost 

Cost 

Rate 

Rate 

F.'i r t No. 

( Days ) 

($) 

( $/unit ) 

( units/day ) 

( units/day ) 

J 


.125 

15 

0.0065 

30000 

.125 

20 

0.1775** 

8000 

.25 

30 

0.1275 

9500 

.125 

10 

0.1000 

7500 

.5 

110 

2.7850 

2000 

.25 

50 

0.2675 

6000 

1 .0 

310 

1 . 5000 

2400 

.5 

130 

5.9000 

1300 

.75 

200 

0.9000 

2000 

.125 

5 

0.0400 

15000 

The  inventory  carrying 

; cost 

is  .10  ($/$/year).  A 

year  consists  of  240 

< 1 : i y : ; . 3 ho  u r s /day . 

In  the  original  problem  [3]» 

this  cost  was  0.1175. 

All  other  authors 

used  0.1775. 

Table  2 

Lotsise  Schedule  for  Doll,  Whybark 
Sequence 


Pi"  nine L 

Production 
Start  Time  (Day) 

Idle  Time 

Following  Production  (Days) 

Beginning 

Inventory  ( Pieces ) 

.125 

0.59 

0.0 

H 

6.21 

0.0 

0.0 

( « 

13.09 

0.0 

0.0 

21.55 

0.0 

0.0 

4 

23.55 

0.59 

0.0 

o 

29.63 

0.0 

0.0 

2 

35.39 

0.0 

0.0 

3 

38.48 

0.43 

0.0 

1 

42.98 

0.0 

0.0 

10 

45.60 

0.0 

0.0 

4 

46.97 

0.59 

0.0 

53.06 

0.0 

0.0 

59.94 

0.0 

0.0 

68.40 

0.0 

0.0 

70.40 

0.59 

0.0 

(\> 

76.48 

0.0 

0.G 

82.73 

0.0 

0.0 

3 

85.33 

0.0 

0.0 

f, 

39.52 

1.55 

0.0 

if) 

92.45 

0.0 

0.0 

4 

93.82 

0.59 

0.0 

99.91 

0.0 

0.0 

Q 

106.78 

0.0 

0.0 

/ 

115.25 

0.0 

0.0 

4 

117.25 

0.59 

0.0 

p. 

123.33 

0.0 

0.0 

129.58 

0.0 

0.0 

132.18 

0.18 

0.0 

7 

1 37. 30 

0.0 

0.0 

10 

1 39. 30 

0.0 

0.0 

140.67 

0.59 

0.0 

146.76 

0.0 

0.0 

9 

153.63 

0.0 

0.0 

162.10 

0.0 

0.0 

4 

164.10 

0.59 

0.0 

170.13 

0.0 

0.0 

176.43 

0.0 

0.0 

179.02 

0.0 

0.0 

C 

133.22 

1.55 

0.0 

10 

186.15 

0.0 

0.0 

I,  $ "2. 071 /day 
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Table  3 

Lotsize  Schedule  for  Practical 
Lotsizes 


Product  Lon 

Idle  Time 

Beginning 

. • ■ : i • t • Product 

Start  Time  (Pay) 

Following  Production  (Pays) 

Inventory  ( Pieces ) 

7 

.125 

0.0 

1715.90 

V* 

4.39 

0.0 

364.63 

9 

10.87 

0.0 

0.0 

/ / 

21.20 

0.0 

0.0 

g 

25.96 

2.13 

0.0 

t t t 

33.45 

0.0 

0.0 

l 

36.95 

0.0 

0.0 

' 10 

39.47 

0.0 

0.0 

o 

41.20 

0.0 

0.0 

■ 0 8 

45.96 

0.0 

0.0 

■ i > 

51.44 

1.65 

0.0 

f, 

58.40 

0.0 

0.0 

L ^ /t 

60.12 

0.0 

1715.90 

64.89 

0.0 

364.63 

L r>  0 

70.87 

0.0 

0.0 

( i 4 

81 . 20 

0.0 

0.0 

i '/  >j 

85.96 

2.13 

0.0 

93.45 

0.0 

0.0 

0 

96.95 

0.0 

0.0 

L0 

99.47 

0.0 

0.0 

1 * 4 

101.20 

0.0 

0.0 

< p 

105.96 

0.0 

0.0 

? '4  3 

111.44 

0.18 

0.0 

V» 

116.80 

0.0 

0.0 

4 

120.13 

0.0 

1715.90 

'♦  > 

124.84 

0.0 

364.63 

1 / v 

130.87 

0.0 

0.0 

141.20 

0.0 

0.0 

><  j 

145.96 

2.13 

0.0 

153.45 

0.0 

0.0 

; 1 r. 

156.95 

0.0 

0.0 

>. 1 in 

159.47 

0.0 

0.0 

161.20 

0.0 

0.0 

165.96 

0.0 

0.0 

171.44 

1 .65 

0.0 

178.40 

0.0 

0.0 

180.13 

0.0 

1715.90 

184.89 

0.0 

364.63 

< • o 

190.87 

0.0 

0.0 

» 0 4 

201.20 

0.0 

0.0 

♦ 1 c 

205.96 

2.13 

0.0 

213.45 

0.0 

0.0 

r 

216.95 

0.0 

0.0 

in 

219.47 

0.0 

0.0 

j 

221.20 

0.0 

0.0 

225.96 

0.0 

0.0 

. 7 '4 

231 .44 

0.10 

0.0 

r • 

237.60 

0.0 

0.0 

‘I.  Vi' 
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